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Resumen 


Actualmente se cuenta con métodos y modelos matemáticos para el 
cálculo y la simulación del ciclo hidrológico, que es una etapa fundamental 
para la correcta contabilización de los recursos hídricos en una cuenca. El 
objetivo de este trabajo es el de realizar un nuevo ensamble de modelos 
de simulación hidrológica y de gestión de recursos hídricos, así como 
evaluar su desempeño contrastando los resultados con los datos oficiales 
reportados en el actual Plan Hidrológico para la subcuenca del río Cabe, 
en Galicia, España. Los resultados de las simulaciones hidrológicas 
obtenidos se compararon con dos estudios similares en la región para 
verificar su validez. Finalmente, en la etapa de simulaciones de la gestión 
de recursos hídricos se encontró que los resultados de este estudio 
presentan una mayor aproximación a los datos observados en 
comparación con los resultados reportados en el Plan Hidrológico vigente. 
Por lo tanto, se considera que el ensamble de métodos y modelos 
propuestos en este trabajo es una contribución para obtener resultados 
de calidad que ayuden a mejorar las situaciones de la gestión actual de 
recursos en la región del río Cabe, así como para hacer frente a los 


impactos del cambio climático en los recursos hídricos. 


Palabras clave: balance hídrico, gestión, simulación, río Cabe, Visual- 
Balan, AQUATOOL. 
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Abstract 


Currently, there are mathematical methods and models for the calculation 
and simulation of the hydrological cycle, which is a fundamental stage for 
the correct accounting of water resources in a basin. The objective of this 
work is to carry out a new assembly of hydrological simulation and water 
resources management models, as well as to evaluate their performance 
by contrasting the results with the official data reported in the current 
Hydrological Plan, for the Cabe river sub-basin, in Galicia, Spain. The 
results of the hydrological simulations obtained were compared with two 
similar studies in the region to verify their validity. Finally, in the water 
resources management simulation stage, it was found that the results of 
this study present a greater approximation to the observed data 
compared to the results reported in the current Hydrological Plan. In this 
way, it is considered that the assembly of methods and models proposed 
in this work is a contribution to obtain quality results that help to improve 
the situations of the current management of resources in the Cabe river 
region, as well as to face the impacts of climate change on water 


resources. 


Keywords: Water budget, management, simulation, Cabe River, Visual- 
Balan, AQUATOOL. 
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Introducción 


La definición de “recursos hídricos” se refiere a las cantidades de agua 
que son utilizadas para las actividades de las sociedades humanas. El ser 
humano no es el único usuario de los recursos hídricos, ya que todos los 
ecosistemas, junto a sus elementos bióticos-abióticos, necesitan del agua 
para su funcionamiento e interacciones. La Organización de las Naciones 
Unidas (ONU, 2018) establece que: “El agua es una condición previa para 


la existencia humana y para la sostenibilidad del planeta”. 


Para poder utilizar y gestionar los recursos hídricos sosteniblemente 
es necesario poder cuantificarlos de manera precisa. De acuerdo con 
Diamond (2005), desde que se tiene conocimiento de las civilizaciones, 
como los egipcios, griegos, romanos, mayas y hasta el presente, su 
supervivencia estuvo marcada por la manera en que gestionaron sus 
recursos naturales. Diamond muestra en su investigación cómo la 
incorrecta gestión de los recursos llevó a la extinción o desplazamiento 
de civilizaciones enteras. En la actualidad y con mayor intensidad desde 
lo establecido en el protocolo de Kioto en 1997, se ha intensificado el 
monitoreo y estudio de los recursos hídricos, al grado de que hoy en día 


se cuenta con técnicas matemáticas asistidas por programas 
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especializados que facilitan el análisis-investigación de los recursos 


hídricos desde nivel de cuenca hasta el nivel de ciudades o pueblos. 


Para abordar los desafíos de gestión de recursos hídricos por lo 
general se necesitan respuestas en dos áreas clave: las que aborden 
cuestiones estructurales, incluyendo la generación y adquisición de datos 
e información, infraestructura y operaciones y mantenimiento; y las 
institucionales, que abarcan cuestiones tales como políticas y fijación de 
precios, o conocimiento e información (GWP 8: INBO, 2009). Por otra 
parte, un enfoque propuesto por la Organización de las Naciones Unidas 
de la Alimentación y la Agricultura (FAO, 2019) es el denominado Water 
Accounting and Auditing (WA8A). Una de las perspectivas del Water 
Accounting está fundamentada en la hidrología. Tal perspectiva 
hidrológica está firmemente basada en el entendimiento del proceso físico 
que gobierna los volúmenes y tasas de cambio en el flujo de agua, así 
como de las existencias en diferentes territorios, y/o bajo diferentes 
condiciones climáticas o régimen de manejo. Otro aspecto crítico del 
Water Accounting es que se debe considerar y evaluar tanto la oferta 
como la demanda de los sistemas de suministro de agua. Del lado de la 
oferta es importante conocer lo más preciso posible la disponibilidad en 
espacio y tiempo de lluvia, agua superficial, agua subterránea y fuentes 
de agua poco convencional. Del lado de la demanda es esencial el 
conocimiento de la demanda de los usuarios, cómo son satisfechas, 
patrones de uso consuntivo y no consuntivo, el nivel de servicio que es 
experimentado por los usuarios y sus beneficios (Batchelor, Hoogeveen, 
Faures, € Peiser, 2017). 
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Los modelos de simulación hidrológica son las herramientas 
utilizadas para investigar los procesos hidrológicos, que van desde 
pequeñas captaciones hasta modelos globales. Cada modelo tiene sus 
propias características únicas y sus respectivas aplicaciones. Algunos de 
ellos son completos y emplean la física de los procesos hidrológicos 
subyacentes, y se distribuyen en el espacio y el tiempo. También, cada 
modelo tiene varios inconvenientes, como la falta de facilidad de uso, los 
grandes requisitos de datos, la ausencia de declaraciones claras de sus 
limitaciones, etcétera (Gayathri, Ganasri, 8 Dwarakish, 2015). 
Sorooshian, Sharma y Wheater (2008) clasificaron los modelos como los 
basados en eventos y los continuos. Los primeros producen información 
solo para periodos de tiempo especificos, mientras que los últimos 
producen información continua. Una de las clasificaciones más 
importantes es el modelo empírico, los modelos conceptuales y los 


modelos de base física. 


Según Islam (2011), en resumen, los modelos para la gestión de 
recursos hídricos se pueden clasificar como: 1) modelos de optimización 
basados en maximización de beneficios económicos; 2) modelos de 
optimización basados en minimización de costos de flujo, y 3) modelos de 
gestión de recursos hídricos basados en sistemas dinámicos. Mayer y 
Muñoz-Hernández (2009) afirman que la simulación y optimización son 
los dos enfoques principales para el modelado de cuencas hidrográficas. 
En la simulación, el comportamiento de los recursos hídricos se simula en 
función de un conjunto de reglas que rigen las asignaciones de agua y la 


operación de la infraestructura; mientras que en la optimización, las 
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asignaciones se optimizan en función de una función objetiva y 


restricciones asociadas. 


La gestión integrada de un sistema de recursos hídricos 
socioeconómico y medioambiental ha suscitado una atención 
considerable, y ha demostrado ser más adecuada para la gestión 
sostenible de los recursos hídricos. Los programas informáticos diseñados 
para el Soporte de Decisión (Sistemas de Soporte de Decisión, SSD), 
involucran los recursos hídricos totales y el consumo de agua, con 
contenido complejo y amplia cobertura. Hoy día, los SSD que ya se han 
desarrollado incluyen principalmente técnicas de optimización, la técnica 
de resolución de conflictos y sistemas adaptativos complejos (Fan, Xu, 
Chen, Li, € Tian, 2020). 


Las herramientas de SSD proporcionan soluciones operativas para 
ayudar a los responsables de la formulación de políticas a abordar los 
complejos problemas ambientales del socio-ecosistema a diversas 
escalas. Además de las funciones centrales que apuntan al análisis de 
decisiones, por lo general incluyen capacidades para modelar y, en 
algunos casos, también para la gestión de procesos participativos. Estas 
herramientas de SSD pueden proporcionar el marco operativo para la 
integración de metodologías y enfoques de modelado de simulación, 


planeación participativa y análisis de decisión (Giupponi 8 Sgobbi, 2013). 


En la zona noroeste de España (Galicia) existen algunos estudios en 
cuencas similares que sirven de referencia para el presente trabajo. En 
específico, se habla de los trabajos de Samper et al. (2005); Raposo, 


Molinero y Dafonte (2012) en Galicia Costa, y Samper y Pisani (2013) en 
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la subcuenca del río Valiñas. La diferencia encontrada en los resultados 
se debe posiblemente a que existe una distinta composición geológica 
entre las subcuencas de Vilasouto y Ponte Pena, en donde casi la mitad 
(47 %) pertenece a grupo de "cuarcitas, pizarras, areniscas y limotitas”, 
a diferencia de la subcuenca de Vilasouto, donde este mismo grupo tiene 
solamente un 13.79 %, que resalta la aparición de calizas, ya que esto 


modifica los parámetros hidrogeológicos del acuífero. 


El objetivo de este trabajo fue plantear una metodología a través 
del ensamble de los modelos matemáticos e informáticos Visual-Balan y 
AQUATOOL, para estimar el balance hídrico en la subcuenca alta del río 
Cabe y evaluar las garantías de abastecimiento en sus demandas de agua. 
A pesar de que en la actualidad existe una planificación hidrológica en 
esta área de estudio, los reportes oficiales aún muestran un sesgo con la 
realidad de los flujos de agua al sobrestimar las cantidades de recursos 


hídricos. 


Materiales y métodos 


Área de estudio 
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Al noroeste de España, en la Comunidad Autónoma de Galicia (CAG), 
dentro la municipalidad de Lugo, se encuentra el Sistema de Explotación 
Cabe (SEC), que se refiere al aprovechamiento resultante de la subcuenca 
del río Cabe, la cual forma parte de la Demarcación Hidrográfica del Miño- 
Sil (DHMS). La superficie de la cuenca del SEC es de 735.18 km?. Durante 
su recorrido de 54.94 km recibe afluentes, siendo importantes el Mao, con 
49.47 km de longitud; Cinsa y Carabelos, con 15.69 y 11.35 km de 
longitud, respectivamente. En todo el SEC hay 200.65 km de longitud de 
masas de agua y 925.05 km de longitud de cauces (MAGRAMA, 2015). 
Uno de los elementos importantes para la gestión de los recursos hídricos 
dentro del SEC es el embalse de Vilasouto, que cuenta con una capacidad 
de 21 Hm?, ubicado en la parte alta de la subcuenca. En la Figura 1 se 
muestra la ubicación del SEC dentro de la DHMS y su punto de unión con 


el río Sil. 
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Figura 1. Ubicación de la subcuenca del río Cabe dentro de la 


Demarcación Hidrográfica Miño-Sil. 


El modelo de simulación hidrológica utilizado en el actual Plan 
Hidrológico (PH) es el modelo conceptual y casi distribuido SIMPA 
(Sistema Integrado para la Modelación del proceso Precipitación 
Aportación), actualizado por el Centro de Estudios Hidrográficos del 
Centro de Estudios y Experimentación de Obras Públicas (Cedex), de 
España. En la Tabla 1 se muestran los datos estadísticos básicos de las 
variables hidrológicas según lo reportado por la CHMS en el PH vigente. 
La precipitación media anual es de 1 041.59 mm y una aportación anual 
media de 386.42 hm3/año, mientras que otras variables hidrológicas para 
el SEC son  evapotranspiración potencial de 604.1 mm/año; 


evapotranspiración real de 518.33 mm/año; recarga de acuíferos de 
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158.66 hm*%/año; caudal superficial de 226.60 hm%3/año y caudal 
subterráneo de 159.82 hm3/año; caudal total de 386.42 hm/año. 


Tabla 1. Estadísticos básicos de precipitación y aportación medios para 
el SEC, serie 1980/81-2011/12. Fuente: CHMS (2015). 


Aportación 
Precipitación Año Año 
Datos anuales Ml ] anual . 
(mm/año) ocurrencia ocurrencia 
(hm*/año) 
Valor mínimo 638.9 2011/2012 140.22 2011/2012 
Valor medio 1041.59 386.42 
Valor máximo 1751.33 2000/2001 841.86 2000/2001 
Coef. var. 0.23 0.40 
Coef. sesgo 0.81 0.86 
Autocorrelación -0.19 -0.22 


En la información reportada en el PH, los usos del agua se dividen 
por “Unidades de Demanda” de tipo urbana (incluye doméstico, industrial, 
municipal, comercial y turismo), agrícola, ganadera y recreacional. En la 
Tabla 2 se pueden observar las cantidades de agua que se extraen de las 
diversas fuentes para satisfacer los 17.78 hm* anuales de todas las 
demandas en el SEC estimadas al 2012. La mayoría del agua se obtiene 
de fuentes superficiales debido a la disponibilidad natural de agua en la 


región. 
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Tabla 2. Origen del abastecimiento en el SEC. Adaptado de CHMS 
(2015). 


Origen de agua en el sistema de explotación 


Sistema de 
Superficial Subterránea Manantial Vol. total 


(hm3/año) (hm3/año) (hm3/año) (hm3/año) 


explotación 


Subcuenca alta del río Cabe 


El área de estudio específica para este trabajo de investigación es en 
realidad una simplificación o reducción del SEC (de 735.18 km), la cual 
se trata de la Subcuenca Alta del Río Cabe (SARC), con un área de 426.34 
km2?, que es aproximadamente la mitad del SEC. El punto de la 
delimitación (salida) de la SARC está dentro de la ciudad de Monforte de 
Lemos, Lugo. El área de estudio se ha establecido a este punto debido a 
que la última estación de aforo (EA1765) se encuentra a 4 km aguas 
arriba en la periferia de la ciudad de Monforte de Lemos. En consecuencia, 
no es posible calibrar un modelo hidrológico después de la estación de 
aforo y realizar una simulación hidrológica en áreas sin calibrar generaría 
resultados con una alta incertidumbre que no son posibles contrastar con 
alguna realidad. Con este fundamento y por tratarse de Monforte de 


Lemos, una ciudad con relevancia en la zona, se estableció el punto de 
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salida de la SARC en las coordenadas latitud 42% 31' 16.23" N y longitud 
7% 30' 56.48" O. La demanda de recursos hídricos por cada tipo en la 
SARC se muestra en la Figura 2, información que ha sido extraída y 
adaptada del PH vigente. Se puede observar que la unidad de demanda 
urbana (UDU) es de 4.2 hm, que representa el 24 % del total y se 
mantiene relativamente estable a lo largo del año. Lo mismo se observa 
para las demandas ganaderas (UDG), demanda recreacional (UDP 10), 
que representan el mínimo de la demanda, al ser del 3 % cada una; la 
unidad de demanda agrícola (UDA), que representa el 70 %, con 12.1 
hm?, y tiene el mayor consumo en los meses de mayo a septiembre y un 


consumo nulo de noviembre a febrero. 
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Figura 2. Usos de los recursos hídricos por cada unidad de demanda en 
la SARC. Adaptado de CHMS (2015). 


Metodología 


La metodología consiste en un ensamble de métodos matemáticos y de 
software especializados para la simulación de procesos del ciclo 
hidrológico y gestión de recursos hídricos. En la Figura 3 se muestra el 


inicio del proceso y el orden en que el ensamble de modelos se va 
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realizando para llegar a los resultados parciales y final del propósito de 
este trabajo, donde cada fase está diferenciada en distinto tono de azul; 
los pasos enmarcados en rojo representan la información de entrada 
obtenida de diversas fuentes consideradas como oficiales; los enmarcados 
en negro hacen mención a software; en amarillo, a modelos matemáticos, 
y los recuadros verdes representan resultados que se consideran 
relevantes para el propósito de este trabajo de investigación. El proceso 
de la Figura 3 inicia con la definición de las subcuencas que componen el 
área de estudio, continuando con la obtención de las series completas de 
los parámetros meteorológicos requeridos por el software de modelación 
hidrológica; después trata de la simulación hidrológica con el software 
Visual-Balan, el cual se alimenta de la información de las dos fases 
anteriores, además de los todos parámetros físicos de la cuenca y de los 
datos de las estaciones de aforo. Así, Visual-Balan realiza una primera 
validación de la simulación hidrológica que, en caso de no cumplir los 
criterios, necesitará una recalibración de los parámetros físicos de la 
cuenca hasta obtener una validación aceptable. Para evaluar el 
comportamiento de modelos, Moriasi et al. (2007) recomiendan el 
número de eficiencia de Nash y Sutcliffe (2003), el coeficiente de 
determinación (R?) y la raíz media del error cuadrático (RMSE). Como 
resultado se obtienen las series de datos de aportación de recursos 
hídricos para cada subcuenca en el área de estudio. La parte que 
corresponde a las simulaciones con el software AQUATOOL requiere de la 
información de los flujos de agua dentro del sistema, como son las 
demandas, retornos, bombeos y red de distribución, entre otros aspectos, 


para posteriormente generar la red de flujos de agua para el escenario en 
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estudio. De acuerdo con Andreu, Solera, Capilla y Ferrer (2007), el 
módulo de SIMGES de AQUATOOL se alimenta de la red de flujos y de las 
series de datos de aportación, realizando la simulación de la gestión para 
abastecer las demandas dentro del sistema. Los elementos de 
almacenamiento (como embalses y reservorios) dentro del sistema de 
flujos sirven de puntos de control para valorar el correcto funcionamiento 
del modelo, a fin de poder comparar los niveles simulados en el embalse 
con los medidos. Como resultado final se obtiene un panorama de la 
contabilización de los recursos hídricos desde antes de su utilización hasta 
la distribución y evaluación del abastecimiento de las demandas para el 


periodo de tiempo evaluado. 


M Información OFICIAL de 
Di en A I Datos observados flujos de agua: demandas, 
pesicin Hsalos retornos, bombeos, 


Terreno e 
canales, caudal ecológico, 


AQUATOOL 


NE .. 


Delimitación Modelación Series de 
subcuencas hidrológica datos de 


VISUAL- aportación Red de bid de agua 
BALAN para el escenario 


Series de datos 
completas | Re- 
T Calibración Simulación de gestión 
í Saca SIMGES-AQUATOOL 
Interpolación - 
espacial Parámetros 
p físicos Subcuenca 
A — Cuantificación recursos 
Parámetros PTE z 
E hídricos, para el periodo 
meteorológicos 
| evaluado 
observados) 
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Figura 3. Proceso de ensamble de la metodología. 


Para la definición de las subcuencas en el área de estudio se requirió 
información geoespacial del Instituto Geográfico Nacional (IGN) de 
España y del software QGIS para procesar y obtener los mapas de 
delimitación de las subcuencas en el área de estudio. Con la información 
de los modelos digitales del terreno (MDT) se trabajaron con el software 
de información geográfica QGIS y las herramientas de GRASS. A través 
del IGN se obtuvieron del centro de descargas la información de los mapas 
con Sistema de Coordenadas de Referencia (SCR) ETRS89, de MDTO5 
LIDAR (Modelo Digital del Terreno y Laser Imaging Detection and 
Ranging) y BTN25 (Base Topográfica Nacional). 


Dentro de los complementos de QGIS que se usaron en este trabajo 
están las utilidades y herramientas de GDAL (Geospatial Data Abstraction 
Library) y GRASS-PROJECT (Geographic Resource Analysis Support 
System). 


En lo que se refiere a las series de datos meteorológicos diarios de 
temperatura media (T media), precipitación (Pr) y evapotranspiración de 
referencia (ETo) para la SARC, se obtienen según los resultados del 
trabajo de Rangel-Parra (2019), donde se detallan los cálculos de la 
estimación por interpolación espacial para las variables climáticas 
mencionadas. 


Los parámetros físicos del suelo del área de estudio se obtuvieron 
en los mapas interactivos de propiedades del suelo de Galicia del 
Departamento de Edafología y Química agrícola de la Universidad de 
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Santiago de Compostela (USC, 2018). Los parámetros físicos que no han 
sido posibles obtener de la bibliografía se ajustaron mediante la 


herramienta de calibración automática que ofrece el Visual-Balan. 


Los datos de aforos observados en los ríos Mao (principal tributario) 
y Cabe se obtuvieron del anuario del Cedex (Cedex, 2018) para las 
estaciones de Ribas Altas, Ponte Pena y embalse de Vilasouto. Para la 
estación de Ponte Pena corresponden al periodo completo de 2007 a 
2017; mientras que para la estación en Ribas Altas solo se utilizaron los 
caudales de invierno de 2006-2017 (cuando no se realiza la actividad de 
regadío), pues durante el verano se hacen extracciones para regadío que 
no están debidamente registradas, por lo que sería posible aplicar la 
restitución al régimen natural. Por último, para el caso del embalse de 
Vilasouto, sí fue posible y necesario aplicar una restitución al régimen 


natural (RRN), al contar con los datos diarios del nivel del embalse. 


De acuerdo con Lopez-Garcia (1993), para calcular el caudal con 
restitución al régimen natural, la ecuación general tiene el siguiente 


planteamiento: 
N=S+T+B-=V-Q=R+EFA (1) 


En donde N es el caudal de un río en régimen natural para el 
periodo; S, el caudal aforado; T, el caudal derivado o trasvasado por 
tomas aguas arriba (hacia fuera de la cuenca); B, el caudal detraído por 


bombeo en acuíferos conectados; V, el caudal vertido por usos aguas 
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arriba; Q, el caudal trasvasado de otras cuencas (hacia dentro la cuenca); 
R, el caudal aportado por recarga adicional en acuíferos; E, el caudal 
evaporado en embalses y masas de agua nuevas, y A es el caudal 


almacenado o liberado por embalses. 


Para este trabajo en particular y aplicado al embalse de Vilasouto, 
se ha modificado la ecuación para adaptarla al caso en particular, 
quedando de la siguiente forma. El caudal total con RRN en el punto de 
salida del embalse se calcula con la siguiente ecuación de balance hídrico 


para el almacenamiento: 

A; =Aj¡-1 +P+N— Evap.—OQsatida (2) 
O para la caudal total: 

N = (A; — Aj¡-1) + Evap. +Osatida — P (3) 


En donde para el día ¡, A; es el almacenamiento en el embalse para 
cierto día y A¡-4 es el almacenamiento del día anterior en unidades de 
volumen/tiempo que corresponde al día ¡; Evap., la evaporación en la 
superficie del embalse; Osalida, el caudal de salida del embalse, y P es la 
precipitación que cae en el área del embalse (se multiplica por el área del 
embalse). La evaporación (Evap.) del embalse de Vilasouto se estimó 
como una fracción de la evapotranspiración de referencia (ETo) calculada 


por interpolación espacial para el centro del embalse. La ecuación 
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reportada por la Asociación Americana de Ingenieros Civiles (ASCE, 1996) 


es: 
ET. = K¿ET (4) 


Hay que agregar el área del embalse para tener el volumen de toda 


la evaporación en el embalse para el día correspondiente: 
K¿(ET, * A) = Evap. (5) 


En donde el valor de Kc es de 0.75 según reportado por ASCE para 
el agua en superficie abierta. Se considera que ET. correspondería a la 
evaporación de la superficie de agua (ET. = Evap.); A es el área del 


embalse, y ETo es la evapotranspiración de referencia en el embalse. 


El software Visual-Balan realiza balances de agua en el suelo edáfico 
en la zona no saturada y en el acuífero, utilizando un entorno interactivo 
de entrada de datos y posproceso de resultados. Los términos 
fundamentales del balance son las entradas por precipitación y riegos; las 
salidas por escorrentía superficial, evapotranspiración, intercepción, flujo 
hipodérmico y flujo subterráneo; la variación del contenido de humedad 
del suelo, de la zona no saturada y del nivel de agua en el acuífero. Estos 
resultados se pueden comparar con datos disponibles de niveles y aforos 
(Samper, 1999). 
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Los criterios que utiliza AQUATOOL para evaluar los fallos en las 
garantías de abastecimiento han sido los de la Instrucción para la 
Planificación Hidrológica (IPH) (IPH, 2008) del gobierno de España para 
demandas urbanas y agrícolas (criterio UTAH DWR): “Para la demanda 
urbana, se considera fallo cuando en el déficit en un mes es superior al 
10 % de la demanda mensual. Para la demanda agrícola, se considera 
fallo cuando el déficit en un año es superior al 50 % de la demanda anual. 
(D.A.). Se considera fallo cuando en dos años consecutivos la suma del 
déficit es superior al 75 % de la D.A., y se considera fallo cuando en diez 
años consecutivos la suma de déficit es superior al 100 % de la D.A.” 


(Solera-Solera, Paredes-Arquiola, 8 Andreu-Álvarez, 2015). 


Resultados 


En la Figura 4 se muestra el resultado del mapa generado en QGIS, en 
donde la Subcuenca Alta del Río Cabe (SARC) se representa en bordeado 
rojo con 426.34 km?, así como de las tres subcuencas internas: Ribas 
Altas, 350.79 km2; Ponte Pena, 35.24 km?, y Vilasouto, 49.24 km?. Las 
subcuencas se establecieron en referencia a las estaciones de aforos 
disponibles. La parte sombreada del mapa corresponde a la subcuenca 
baja del río Cabe (308.84 km?, que no se incluye en este trabajo), la cual 
se une al río Sil hacia el sur. Las coordenadas geográficas WGS84 de sus 
centroides son SARC 420 36' 13” N y 7% 25' 2.9” W; Ribas Altas 420 37” 
40” N y 79 23" 40” W; Embalse Vilasouto 42% 41' 10” N y 70 21” 12” W; 
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por último, Ponte Pena 42% 39" 10” N y 7% 18' 37” W, con proyección 
ETRS89/UTM zona 29 (EPSG:25829). 


Embalse de Vilasouto 


? 
) ) 
( A 
| 


Leyenda 


"| 0 Estaciones de aforo 

E Subcuenca Ponte Pena 
EJ Subcuenca Vilasouto 

E Subcuenca Ribas Altas 
(A Subcuenca Alta río Cabe 
Subcuenca del río Cabe 
DH-MiñoSil 


Figura 4. Resultado de la generación de subcuencas dentro de la 


Subcuenca del río Cabe dentro. 


Respecto al cálculo de la restitución al régimen natural (RRN), en la 
Figura 5 se puede observar el flujo promedio diario de salida en el embalse 


de Vilasouto restituido. Al graficar este caudal con RRN, se observó que 
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después del día 20/06/2013 de los 11 años hidrológicos había un 
comportamiento irregular negativo en los caudales (círculos indicados) 
debido al error en la estimación de los caudales de salida, y en sus 


registros y validación. Por tal motivo, solo se utilizaron los datos 


anteriores a esta fecha de la irregularidad, que corresponde del 
01/10/2006 al 19/06/2013, como referencia para calibrar el modelo 


hidrológico Visual-Balan (puntos azules de la Figura 6) en esta subcuenca. 


m3/s 


Figura 5. Resultado de la RRN en el embalse de Vilasouto, flujo 


promedio por día. 
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—— Caudales calculados (m3/s) 
» Caudales medidos (m3/s) 


—— Precipitación (mm) 


CAUDALES PRECIPITACIÓN 
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Figura 6. Comparación de la calibración en Vilasouto con los caudales 


18/01/10 
24/02/11 
01/04/12 
07/05/13 
20/07/15 
24/08/16 
30/09/17 


observados y simulados. 


En la Figura 6 se muestra el resultado de la calibración en la 


subcuenca de Vilasouto, comparando los caudales observados y 
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simulados después de haber eliminado el periodo con flujos negativos 


explicado en la Figura 5. 


En la Tabla 3 se muestra el resultado de aportaciones media anual 
por subcuenca en la SARC para el periodo 2006-2017; la subcuenca de 
Vilasouto, 29.71 hm*/año (12 %); Ponte Pena, 24.26 hm*%/año (8 %); 
Ribas Altas, 132.32 hm3/año (62 %); el área restante entre Ribas Altas y 
el punto de salida de la SARC aporta 29.88 hm*/año (18 %), para un gran 
total de 216.36 hm*%/año. Según lo reportado por la CHMS, el SEC 
completo aporta 386 hm“/año, por lo que los 216.36 estimados en este 


trabajo para el SARC representan el 56 % de lo estimado para el SEC. 


Tabla 3. Aportación de recursos hídricos por subcuenca de la SARC. 


Aportación 
(hm3/año) 


Ribas Altas 132.32 
Remanente SARC 29.88 
Total en la SARC 216.36 100 


Después de estimar el caudal con RRN, se hizo la simulación 


hidrológica con Visual-Balan en la subcuenca de Vilasouto de forma 


iterativa y recalibrando los parámetros físicos del suelo, utilizando las 
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ayudas de las herramientas incluidas en Visual-Balan para análisis de 
sensibilidad y autocalibración. Los parámetros físicos del suelo del área 
de estudio se obtuvieron en los mapas interactivos de propiedades del 
suelo de Galicia del Departamento de Edafología y Química Agrícola de la 
Universidad de Santiago de Compostela (USC, 2018). En la Figura 6 se 
muestran los caudales simulados con el modelo calibrado y los caudales 
observados de la base de datos del anuario del Cedex, en donde se obtuvo 
un valor del coeficiente de Nash-Sutcliffe-Efficiency (NSE) de 0.69. Los 
caudales simulados se ajustan adecuadamente a la magnitud y tendencias 
de los datos observados. Por otra parte, el resultado de la simulación 
hidrológica para Ribas Altas obtuvo una validación de NS de 0.82, y para 
la subcuenca de Ponte Pena se tuvo un valor de NS de 0.38. Moriasi et al. 
(2007), en su publicación, dieron una clasificación de interpretación a los 


valores del NSE, los cuales se muestran en la Tabla 4. 


Tabla 4. Interpretación de los valores del criterio Nash-Sutcliffe- 


Efficiency. 
Calificación del rendimiento NSE 
Muy buena 0.75 < NSE > 1 
Buena 0.65 < NSE > 0.75 
Satisfactoria 0.5 < NSE > 0.65 
Insatisfactoria NSE =< 0.5 


Fuente: Moriasi et al. (2007). 
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El NSE obtenido en la subcuenca de Ponte Pena es el más bajo de 
las tres subcuencas simuladas, no obstante, solo representa el 8 % en las 
aportaciones del área de estudio, como se muestra en la Tabla 3. Por otra 
parte, no se encontró información detallada de los parámetros geológicos 
y de flujos subterráneos para mejorar la validación de la simulación en 
Ponte Pena. Para aumentar la consistencia de la calibración obtenida para 
la SARC (y observar el impacto de resultado obtenido en Ponte Pena), se 
contrastó la distribución de flujos con otros dos estudios hidrológicos 
dentro de la zona de Galicia (Raposo et a/., 2012) para la Demarcación 
Hidrográfica, y para el río Valiñas (Samper 8 Pisani, 2013). En la Tabla 5 
se puede observar la similitud en la importancia de los componentes 
hidrológicos. Además, en la Figura 7 se muestra la distribución gráfica 
entre los dos resultados de flujos del caso del río Valiñas y la SARC. La 
coincidencia de tales comparaciones respalda los resultados obtenidos 
para el área de estudio de este trabajo y la representatividad en la 


calibración obtenida con Visual-Balan. 
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Tabla 5. Comparación de los resultados de Visual-Balan para la SARC 


con dos estudios en la misma zona de Galicia. 


A ] EL Río Valiñas Galicia Costa 
Parámetros hidrológicos ] ] 
Samper y Pisani Raposo el al. 
(valores anuales) 
(2013) (2012) 


Flujo subterráneo (mm) 


Caudal total (mm) 
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Escorrentía Flujo Flujo 
superficial ETR (mm) hipodérmico | subterráneo 
(mm) | (mm) | (mm) 


Precipitación | Interceptación 
(mm) (mm) 


Caudal total 
(mm) 


C Río Valiñas 


100.00% 15.08% 0,46% 34.78% 40.43% 9.07% 


50.05% 


| mCabe_R 


100.00% 9.68% 1.69% 39.33% 39.84% 9.49% 


51.03% 


Figura 7. Comparación de los resultados de las variables hidrológicas la 


SARC y la subcuenca del río Valiñas. 


Una de las características de Visual-Balan es la capacidad de realizar 
una doble calibración frente a caudales superficiales y niveles 
piezométricos. Sin embargo, en lo que respecta a los recursos hídricos 
subterráneos, no se cuenta con información sobre los datos de niveles 
piezométricos en el área, información que es fundamental para mejorar 
los resultados de calibración del modelo hidrológico-geológico. De nuevo 
Raposo et al. (2012), y Samper y Pisani (2013) coinciden en sus trabajos 
con el reto de calibrar el modelo hidrológico al trabajar en un área de 
estudio con variaciones en la geología de la cuenca, y la heterogeneidad 
de las características de almacenamiento y transmisión del agua 


subterránea en rocas de tipo cristalino. 
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Para la simulación de gestión en AQUATOOL, las aportaciones 
calculadas en Visual-Balan para cada subcuenca se consideran como la 
información de entrada de recursos hídricos. El esquema de flujos 
desarrollado en AQUATOOL para la SARC de la Figura 8 es la 
correspondiente adaptación del existente en el PH 2016-2021 para el SEC. 
Adicionalmente fue agregado el caudal ecológico (establecido igualmente 
en el PH para esta subcuenca) en forma de demanda de recurso hídrico; 


así fue posible evaluar su garantía de cumplimiento. 
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El nivel de reserva observado en el embalse de Vilasouto fue 
considerado el punto de referencia para evaluar la simulación de flujos de 
AQUATOOL contra el comportamiento real del sistema. En la Figura 9 se 
muestran la gráficas entre los niveles observados y simulados por 
AQUATOOL para la SARC, así como los niveles de la simulación llevados 
a Cabo y reportados por la CHMS en el PH vigente. Se puede observar una 


clara mejor aproximación a la realidad de los resultados de este trabajo 


que los reportados en el PH. 
j | 
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Figura 9. Comparación entre los niveles del embalse de Vilasouto 


simulados con AQUATOOL y niveles observados. 
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La diferencia encontrada entre los resultados de la CHMS y los de 
este trabajo recae en las estimaciones de aportación de recursos hídricos. 
Mientras que los datos utilizados por la CHMS fueron generados por el 
software SIMPA, en este trabajo se utilizó la información de los 
parámetros meteorológicos reportados en Rangel-Parra (2019) derivados 
de los datos de MeteoGalicia. La simulación hidrológica en Visual-Balan 
implica un proceso de mayor complejidad que el modelo SIMPA, por lo 
que las aportaciones generadas por SIMPA están sobrestimadas y al ser 
simuladas en AQUATOOL resultan en una ausencia de fallos en las 


garantías. 


Por último, para los fallos en la garantía de abastecimiento se 
utilizaron los criterios del IPH preestablecidos en AQUATOOL. En la Tabla 
6 se muestra el resumen de los resultados de los fallos para el periodo 
simulado 2006-2017. Se encontró un total de 10 fallos para las 
demandas; el caudal ecológico es el de mayor número con tres; la 
demanda ganadera en Bóveda presenta dos, y por último las que 
presentan un fallo fueron la demanda agrícola del valle de Lemos, la 
demanda ganadera en Monforte, la demanda agrícola en Monforte, la 
demanda agrícola en Bóveda y la demanda recreacional UDP10. Para las 
urbanas (UDU) no se presenta ningún fallo. Es importante señalar que a 
las UDU se les ha dado prioridad sobre el resto de las demandas. A pesar 
de los fallos encontrados por los criterios del IPH, la garantía en 
porcentaje volumétrico se mantiene arriba de 99 %. Las fechas de los 
fallos coinciden con los periodos de menor nivel de reserva en el embalse 


de Vilasouto. 
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Tabla 6. Fallos en la garantía de abastecimiento para las demandas de 
la SARC en la simulación con AQUATOOL. 


Núm. , , z Max. déficit, 
Garantía | Garantía | Max. déficit, un 
Demanda de dos meses 
(%) Vol. (%) mes (hm) 
fallos (hm) 
0.034 0.034 
0.021 0.023 


UDA_3301_Valle de 


Lemos 


UDG Ficticia 2301 y 
2302 Bóveda 


UDG ficticia 2303 


Monforte de Lemos 


ML 
UDA ficticia 2303 
1 99.2 0.011 0.011 
Monforte de Lemos 


UDA Ficticia 2301 y 
2302 Bóveda 


0.045 0.045 


0.032 0.032 


UDP 10 


Caudal Ecológico 
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Conclusiones 


La aplicación conjunta de los resultados de caudales obtenidos con el 
modelo hidrológico Visual-Balan y del modelo de soporte a la decisión 
AQUATOOL produce valores más próximos a los valores registrados en el 
embalse de Vilasouto que los reportados en el Plan Hidrológico vigente 
para el río Cabe. Por lo tanto, la mejora en la simulación del ciclo 
hidrológico de este trabajo se considera una contribución para los futuros 


planes hidrológicos en la Demarcación Hidrográfica del Miño-Sil. 


La coincidencia en magnitud y distribución de los flujos de los 
componentes hidrológicos con otros dos estudios en la misma región 
respalda los resultados de la calibración del modelo hidrológico Visual- 


Balan. 


En la simulación con AQUATOOL se encontró que existen fallos en 
las garantías de abastecimiento para algunas demandas de agua que no 
están identificadas en el Plan Hidrológico. Estos fallos coindicen con los 
niveles bajos en el embalse de Vilasouto. El reporte en el Plan Hidrológico 
se ve afectado por la sobrevaloración de los recursos hídricos del modelo 
SIMPA, que se reflejan en el embalse con niveles mayores a los reales, 
resultando en menor estimación de fallos en el abastecimiento de las 


demandas de agua. 


La ausencia de datos de niveles piezométricos es una limitación en 


el momento de calibración de modelos hidrológicos, sobre todo en una 
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zona como la de Galicia y su naturaleza de hidrogeológica. Sin embargo, 
con los datos disponibles actualmente es posible lograr una calibración 


funcional para simular los flujos de agua superficiales. 


Contar con este ensamble de modelos Visual-Balan y AQUATOOL 
calibrado y validado abre la posibilidad de extender las investigaciones en 
la subcuenca alta del río Cabe para evaluar los posibles efectos del cambio 
climático sobre los recursos hídricos. Con la ayuda del ensamble realizado 
en este trabajo se puede abordar con mayor énfasis cualquier tema en la 
gestión de los recursos hídricos y el abastecimiento para las demandas de 


agua. 
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